Excited Thomas-Efimov levels in ultracold gases 
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Since the early days of quantum physics, the complex behavior of three interacting particles has been the 
subject of numerous experimental and theoretical studies. In a recent Letter to Nature, Kraemer et al. [Nature 
(London) 440, 315 (2006)] report on experimental "evidence for Efimov quantum states" in an ultracold gas of 
cesium atoms. Such quantum states refer to an infinite series of energy levels of three identical Bose particles, 
accumulating at the threshold for dissociation as the scattering length of each pair is tuned to infinity. Whereas 
the existence of a single Efimov state has been predicted for three helium atoms, earlier experimental studies 
concluded that this elusive state had not been found. In this paper we show by an intuitive argument and full nu- 
merical calculations that the helium and cesium experiments actually provide evidence of the same, ground state 
of this trimer spectrum, which the helium experimentalists and pioneering theoretical studies had not associated 
with Efimov's effect. Unlike the helium trimer, the observed l33 Cs 3 resonance refers to a Borromean molecular 
state. We discuss how as yet unobserved, excited Efimov quantum states might be detected in ultracold gases of 
85 Rb and of 133 Cs at magnetic field strengths in the vicinity of 0.08 T (800 G). 

PACS numbers: 34.50.-s,03.75.-b,34.10.+x,21.45.+v 



I. INTRODUCTION 

Efimov's effect HI! refers to a scenario in which three 
identical Bose particles interact via weak pair potentials, each 
supporting only a single bound state when the s-wave scatter- 
ing length a is positive. None of the pairs is bound when a 
is negative, while at the intersection of the two regimes the 
scattering length has a singularity. As early as 1935 Thomas 
(4J] showed that under these conditions three particles can be 
tightly bound, despite their weak binary interactions. Thirty 
five years later, Efimov (]]] discovered a striking extension of 
Thomas' effect, predicting the existence of an infinite series 
of excited three-body energy levels in the limit \a\ — > oo. Each 
Efimov state emerges when a is negative at the threshold for 
dissociation into three free particles and vanishes eventually 
via decay into a two-body bound state and a third free parti- 
cle as a is increased to positive values across the singularity. 
Such an ideal universal scaling behavior of three-particle en- 
ergy levels is illustrated in Fig.Q] Two of the levels (labeled 
by their degree of excitation, n = 1, 2, 3, . . .) are visible but 
ideally they all have the same functional form as a is varied, 
the full extent of which is displayed only for the E„ level. 

Strict universality in a three-body system presupposes a 
large scattering length whose modulus by far exceeds all the 
other length scales set by the pair potentials. This specific 
requirement on the two-particle interactions in combination 
with the wide spread of energy scales in Fig.Q]have for a long 
time prevented any generally accepted discovery of Efimov 
states in nature. Experimental realizations of Efimov's sce- 
nario have been suggested, for instance, for trimer molecules 
consisting of three 4 He atoms [5], as well as for alkali atoms 



in ultracold gases with magnetically tunable interactions (see 
Refs. H, 0, W 

and references therein). Helium trimers well 
match the assumptions underlying Thomas' and Efimov's ef- 
fects, as for each atom pair there exists only a single weakly 
bound dimer state giving rise to the large positive scattering 
length of about [9] 200 oo (flo being the Bohr radius). The 
inter-atomic interactions are inaccessible to external manipu- 
lation but are sufficiently weak to give rise to two trimer bound 
states of the Thomas-Efimov spectrum. 

The excited 4 He3 Ei level follows the qualitative trends in- 
dicated in Fig. Q] when the pair potential is scaled to mimic 
a hypothetical tuning of a. In particular, it vanishes into the 
atom-dimer continuum as the pairwise attraction is increased 
fTol [Till and has therefore been referred to as a genuine Efi- 
mov state. By contrast, a rigorous variational estimate lfl2tl 
reveals that the binding energy of the ground state of the 
Thomas-Efimov spectrum (not shown in Fig. Q} violates the 
global universal scaling behavior of Fig. Q] as it never inter- 
sects with the two-body level. Similar deviations from three- 
body universality, as \Ja is hypothetically increased, have 
been discussed in Ref. [ 13], in the context of the three-nucleon 

JUS 



problem. In particular, several publications JlC 



11, 



had not associated the experimentally observed 1 14lll5ll He^ 
ground state with Efimov's effect. In the limit of zero inter- 
action, i.e. l/a — > -00 in Fig.Q] however, this ground trimer 
level must vanish in a three-body zero-energy resonance. 

In this paper we show that recent experiments in ultra- 
cold gases of 133 Cs atoms |§] also provide evidence for the 
ground trimer level of the Thomas-Efimov spectrum via ob- 
servation of its associated zero-energy resonance. We inves- 
tigate several Efimov spectra of the experimentally relevant 
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FIG. 1: (Color online) Efimov plot Qj] illustrating the wave num- 
ber K = -(m\E n \/h 2 ) 1 ' 2 (n = 1,2,3,...) indicated by solid curves 
vs. the inverse scattering length, 1 /a. Here m is the single particle 
mass and E„ is the energy of the nth three-particle level. The wave 
number labeled dimer (dotted line) is associated with the two-body 
bound state emerging in the limit I /a — > (i.e. \a\ — * oo) at the ver- 
tical dashed line. In accordance with Ref. [J], its energy is universal, 
i.e. E b = -h 2 /(ma 2 ). The limits l/a — » +oo of the horizontal axis 
refer to zeros of the scattering length. Ideally, each Efimov level 
emerges in a three-body zero-energy resonance at negative a and 
eventually vanishes into the continuum leading to an atom-dimer res- 
onance at positive a. The ground state of the Thomas-Efimov spec- 
trum (n = 1) is not displayed. 



We discuss prospects for their detection in comparison with 
a previous suggestion for experiments using 85 Rb [20]. Sec- 
tion [IV] summarizes our conclusions. Finally, the appendix 
provides a detailed description of our general method lf23ll 
for calculating resonance-enhanced three-body recombination 
loss-rate constants including the regimes of negative as well as 
small positive scattering lengths. 



H. RESONANCE-ENHANCED THREE-BODY 
RECOMBINATION IN CESIUM GASES 

Whereas the helium studies of Refs. fTil [l5ll relied upon 
direct observation of 4 He3 in a molecular beam, the existence 
of trimer levels in ultracold gases 18J] has been inferred from 
three-body zero-energy resonances as illustrated in Fig. Q] 
Their signatures manifest themselves, for example, in a res- 
onant enhancement of three-body recombination loss [6, 7]. 
In such processes three initially free atoms collide to form a 
dimer, while the binding energy is transferred to the relative 
motion of the third atom with respect to the bound pair. Their 
final kinetic energies are usually high enough for all three col- 
liding particles to leave an atom trap confining an ultracold 
gas, whose atom number N thus decays in accordance with 
the rate equation 



N = -K 3 (n 2 )N. 



(1) 



Here (n 2 ) is the mean square density and Kj, is the three-body 
recombination loss-rate constant. 



alkali atomic species 133 Cs and 85 Rb in the vicinity of sin- 
gularities of the diatomic scattering length, realized by mag- 
netically tunable Feshbach resonances Im ITtL [l8ll . This in- 
volves predictions of the energies of metastable trimer levels 
as well as numerical calculations of three-body recombination 
loss-rate constants yg, iTJ, motivated by the observations re- 
ported in Ref. H]. Our studies indicate that ultracold gases in 
the presence of broad, entrance-channel dominated diatomic 
zero-energy resonances and negative background-scattering 
lengths 11911 may be best suited for the detection of as yet 
unobserved excited trimer levels of the Thomas-Efimov spec- 
trum JH. 

The paper is organized as follows: In SectionlHlwe briefly 
discuss the technique of detecting three-body zero-energy res- 
onances |0, in ultracold gases of 133 Cs at low magnetic- 
field strengths. Based on quantitative calculations of associ- 
ated three-body recombination loss-rate constants, we inter- 
pret the measurements of Ref. [8] in terms of their relation to 
Efimov's effect [1,2]. In this context, we discuss crucial short- 
comings of purely universal treatments of near-resonant three- 
body recombination into comparatively tightly bound target- 
dimer states at negative as well as small positive scattering 
lengths H HI], |22l] . Using an intuitive argument, we illustrate 
our finding that the experimentally observed three-body zero- 
energy resonance [8] refers to the ground state of the Thomas- 
Efimov spectrum. In Section [HI] we predict magnetic-field 
strengths associated with excited state three-body zero-energy 
resonances in ultracold 133 Cs gases in the vicinity of 800 G. 



A. Ultracold collisions of cesium atoms at low magnetic-field 
strengths 

Observation of atom loss was used in Ref. J^] to track down 
a three-body zero-energy resonance in gases of 133 Cs prepared 
in the (F = 3, nif = +3) Zeeman ground state at temperatures 
> 10 nK. Here F refers to the total atomic angular momen- 
tum quantum number, and nif indicates the associated projec- 
tion with respect to an applied, spatially homogeneous mag- 
netic field of strength B. This experimental setup [8] allows 
magnetic tuning of the s- wave scattering length a to mimic 
the variation of inter-particle interactions in Efimov's scenario 
UIH], taking advantage of the occurrence of a diatomic zero- 
energy resonance in cesium collisions at low fields. 

In the vicinity of any diatomic zero-energy resonance [24], 
the scattering length assumes all values between -00 and +00. 
The associated strong magnetic-field variation a{B) can usu- 
ally be described, to an excellent level of accuracy, by the fol- 
lowing general formula l25ll : 

a(B) = a A l -lTB-} (2) 

Here Bq is the position where a has a singularity, a\, g is 
the background scattering length observed asymptotically far 
from Z?o, while the resonance width, AB, refers to the dis- 
tance between Bo an d the zero-crossing point of the scattering 
length, i.e. a(B + AB) = 0. 
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FIG. 2: (Color online) Comparison between full coupled-channels 
(circles) and two-channel (curves) calculations of the I33 Cs 6i-dimer 
energies (upper panel) (l7Ll25ll27ll . as well as their associated closed- 
channel admixtures 1131 . Z(B) (lower panel), vs. the magnetic-field 
strength, B. For each magnetic-field strength zero energy is chosen 
to coincide with the scattering threshold. The resonance position, 
So, is indicated by the dotted line. Curves labeled E i l and E d 2 refer 
to the two highest excited vibrational bound states for B > B and 
the levels they adiabatically correlate with at B < Bo, whereas the 
Feshbach molecular energy E b existing for B < B is not resolved in 
the upper panel. Its closed-channel admixture in the lower panel is 
negligible. 



In the particular case of 133 Cs atoms prepared in the (F - 
3,nif = +3) state near zero field strength, the strong in- 
dependence of the scattering length is driven by a diatomic 
zero-energy resonance in the mirror-spin-channel of atoms in 
the excited (F = 3,nif = -3) Zeeman state. This scenario 
can be formally described ll26ll by a negative resonance posi- 
tion that we calculate at B = -1 1 .2 G (1 G=10~ 4 T) using the 
complete two-body Hamiltonian model of Ref. [27], which 
accurately represents the positions of all known diatomic reso- 
nances in cesium collisions. An associated vibrational molec- 
ular energy spectrum is illustrated in the upper panel of Fig. [2] 
for isotropic 6s-dimer states ifnl |2 5[ |2 7tl. whose highest ex- 
cited, Feshbach molecular, level E\, [19] (unresolved in Fig.[2]i 
causes the zero-energy resonance at Bo- The large, posi- 
tive background-scattering length of a\, g = 1720 ao. in com- 
bination with the resonance width of AB = 27.5 G, allow 
a(B) to be varied between about -2500 ao and zero within a 
field range between B — and about 17 G. Beyond its zero- 



crossing point the scattering length rises to ab g , as shown in 
the inset of the upper panel of Fig. [3] 



B. Three-body zero-energy resonances 
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FIG. 3: (Color online) Calculated three-body recombination loss-rate 
constants A" 3 (upper panel) and energies of metastable trimer levels 
of 133 Cs atoms labeled E t and Ei (lower panel) vs. the magnetic- 
field strength B. For each magnetic-field strength zero energy co- 
incides with the scattering threshold for three asymptotically free 
atoms. Only ranges of B > are accessible to experiments with 
atoms in the Zeeman ground state, whereas B < would be accessi- 
ble to ex peri ments with the (F = 3,m F = -3) level. Measurements 
of K3 1 8, 28] performed at 10 nK and about 200 nK are indicated by 
circles and squares, respectively. The dashed and dot-dashed curves 
labeled £ d , and in the lower panel refer to the 6s-dimer lev- 
els I17L |25L 12711 of Fig. [2] constituting the target states included in 
our calculations. Dashed and dot-dashed curves in the upper panel 
indicate their associated contributions to the total three-body recom- 
bination loss-rate constant (solid curve). The inset of the upper panel 
illustrates the variation of the scattering length with magnetic field. 

Based on numerically exact three-body calculations with- 
out fit parameters using the method of Ref. i23"ll . Fig.[3]shows 
the B-dependence of ^3 as well as the metastable 133 Cs3 
level causing a three-body zero-energy resonance observed 
in Ref. [8] at 7.2 G. Our estimates of its position rely upon 
two independent methods: The first approach II20I1 follows 
Thomas' and Efimov's original treatments and accounts only 
for the weakly bound Feshbach molecular dimer level 111911 . 
which causes the diatomic zero-energy resonance at Bo in the 
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mirror-spin-channel. This leads to a zero crossing of the en- 
ergy of a trimer state at 7.5 G (solid curve labeled E\ in Fig. [3}. 

Our second approach consists in directly determining the 
peak positions in the calculated magnetic-field dependent 
three-body recombination loss-rate constants M23I1 . In addi- 
tion to the Feshbach molecule, these calculations account for 
the next two highest excited 6s-dimer levels [271 below the 
dimer dissociation limit, including their admixtures from dif- 
ferent atomic Zeeman states [19], as discussed in the ap- 
pendix. The associated pair potentials are based on a para- 
metric two-channel representation of the diatomic low-energy 
spectra resulting solely from coupled-channels calculations. 
The one approximation we make in the two-body coupled- 
channels bound-state calculations is to neglect the relatively 
weak coupling of the 6s-levels to nearby dimer levels of pre- 
dominantly (/-wave symmetry. We plan future calculations 
in which we include such levels. Our exact solutions of the 
three-body Schrodinger equation have been performed at zero 
kinetic energy of the three initially free colliding atoms and 
lead to a resonance peak in K3 at 5.1 G. We note that the 
second approach could also be used to calculate the Thomas- 
Efimov trimer levels of Fig. [3] in terms of resonant enhance- 
ments of atom-dimer collision-rate constants [23], rather than 
the stable bound states of the first approach. In such a more 
precise treatment, the atom-dimer resonances become degen- 
erate with the zero-energy scattering threshold exactly at those 
magnetic-field strengths where three-body zero-energy reso- 
nances occur in the recombination loss-rate constant A/?. 

Relative to the width of the diatomic zero-energy resonance 
the deviations between the calculated and observed peak po- 
sitions are less than 8 % for both methods without any adjust- 
ments to the measured Kj,. We note that reducing the temper- 
ature of the cesium gas tends to shift the peak towards lower 
fields 112911 . The comparison between both methods and the ex- 
periments indicates that the three-body zero-energy resonance 
positions in the upper panel of Fig.[3]are largely insensitive to 
the details of the diatomic molecular target levels. This be- 
havior is physically sensible, as these positions depend only 
on the metastable Thomas-Efimov levels close to the three- 
body dissociation threshold in the lower panel. 



Incompleteness of universal approaches to three-body 
recombination at negative scattering lengths 



As opposed to universal fitting procedures 112 U 12211 . our cal- 
culations also correctly incorporate the two- and three-body 
physics in the vicinity of the minimum of A3 (just below 20 G 
in the upper panel of Fig. [3j and fully recover its measured 
position. While the overall calculated magnitudes are in good 
agreement with experiment throughout, we observe a strong 
sensitivity of £"3 to the B-dependent closed-channel admix- 
tures Ill9ll of the diatomic target states illustrated in the lower 
panel of Fig. [2] This sensitivity at negative as well as small 
positive scattering lengths is to be expected because the for- 
mation of a dimer from initially free atoms should be sup- 
pressed when the atoms have to change their Zeeman states 
during this process. These findings imply that depends on 



B not only through a(B) but also through the strongly varying 
properties of the comparatively tightly bound diatomic molec- 
ular target levels of Fig. [2] Consequently, any fit analysis of 
the measured data in Fig.[3]incorporating the /^-dependence of 
K3 only through the scattering length is incomplete, and could 
therefore lead to unreliable conclusions about three-body uni- 
versality in the system. 



D. Interpretation of the cesium experiments 
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FIG. 4: (Color online) Comparison between Efimov plots associated 
with 133 Cs in the vicinity of zero field (upper panel) and 4 He (lower 
panel). The filled circles refer to the degree of excitation (n = 1) 
of the trimer states detected in the experiments of Refs. 0| [l4Tl at 
the inverse scattering lengths indicated on the horizontal axis. Dot- 
ted curves indicate energies of dimer levels causing diatomic zero- 
energy resonances at 1 ja = 0. The range of scattering lengths in the 
upper panel extends from the observed [8] three-body zero-energy 
resonance position to a = at field strengths B < 0. In accordance 
with Refs. fl2l flUl . the curves labeled E[ illustrate the increasing 
separation of the trimer ground level from the dimer energy as 1 ja 
is increased. By contrast, the E 2 level in the lower panel approaches 
the atom-dimer threshold. It eventually gets overrun by the dimer 
energy causing an atom-dimer resonance as the potential is scaled to 
enhance the pairwise attraction flol [TUl . For this reason, only the 
excited E 2 level has been referred to in Refs. (3, QU [Hi [T^l as a 
genuine Efimov state. 

Figure [3] shows that according to our calculations the reso- 
nance peak experimentally observed at 7.2 G [8] is associated 
with the ground trimer level of the Thomas-Efimov spectrum, 
E\. As opposed to the peak positions and magnitudes of £3, 
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the labels attached to the trimer levels in the lower panel are 
independent of the detailed physics included in the approach. 
All that is needed for this classification is the knowledge that 
the observed resonance QSQ is the first to appear on the low- 
field side of the zero crossing of the scattering length. This is 
immediately evident from the I33 Cs Efimov plot in the upper 
panel of Fig. [4] in which the field strength associated with 
a — corresponds to the limit I /a — > -oo. As is de- 
creased the inverse scattering length l/a(B) increases. Zero 
magnetic-field strength is indicated by the vertical dashed 
line. The trimer energy spectrum beyond this point is indi- 
cated by dashed curves in the upper panel of Fig. |4]and might 
be accessed in gases of I33 Cs atoms prepared in the excited 
{F = 3, niF = -3) Zeeman state. 

For comparison, the lower panel of Fig. |4]shows the Efimov 
plot associated with 4 He3 11141 1 1 311 - Both the cesium [8|] and 
helium lfl4l [l5ll experiments provide evidence for the same 
level of the Thomas-Efimov spectrum. Whereas the 4 He3 
molecule is stable with respect to spontaneous dissociation, 
the I33 Cs3 state can decay into an atom and, e.g., one of the 
dimer levels in the lower panel of Fig. [3] in accordance with 
energy conservation. We note in this context that all stable 
133 Cs3 molecular levels are below the diatomic ground-state 
energy, whose typical modulus for alkali systems exceeds the 
largest hx 100 MHz energy scale considered in this paper by 
six orders of magnitude. Consequently, the ground state of the 
Thomas-Efimov spectrum of cesium does not coincide with 
the stable molecular ground state of I33 Cs3. The resonance 
signature of the observed [8] metastable Thomas-Efimov ce- 
sium trimer, however, has been detected at a negative scat- 
tering length. To our knowledge, this measurement therefore 
provides the first experimental evidence for the existence of 
a Borromean rf30Tl molecular state, i.e. a three-body energy 
level, with a finite width in the present context, in the ab- 
sence of any diatomic levels energetically above it. Among 
the exotic three-atom bound molecular states of helium, in- 
cluding the isotope mixture 4 He2 3 He, for instance, at least one 
atom pair, namely 4 He2, is bound 11411 . Accordingly, only a 
"pseudo-Borromean" four-atom complex, 4 He2 3 He2, has been 
observed in molecular-beam experiments ||3lll . 



HI. EXCITED THOMAS-EFIMOV TRIMER LEVELS 

The prospect of an experimental verification of Efimov's 
effect depends on the criterion for its discovery, since a mea- 
surement of an infinite number of trimer levels with increas- 
ingly loose bonds is impossible. Numerous pioneering stud- 
ies, such as Refs. JH OH EL therefore agreed that the de- 
tection of a single Efimov state would be sufficient. Refer- 
ences 1I22I l32tl suggest that the 4 He3 ground state should be 
referred to as the first example of an Efimov state in nature, 
as the modulus of the diatomic scattering length by far ex- 
ceeds all the other length scales set by the pair potentials [9]. 
Several studies summarized in Refs. O20l 13211 have shown that 
the 4 He3 molecule is indeed universal, i.e. both levels are ac- 
curately determined by the scattering length in addition to a 
single independent parameter 01 HH). Such a parameter could 



be provided, e.g., by recent experimental studies lfl5ll showing 
that the bond lengths of the observed 4 He3 molecules Ill4ll are 
consistent with predictions for the ground state. Several theo- 
retical studies suggest that the ground-state triton 3 H nucleus 
113 311 might also be described using universal treatments (see, 
e.g., Refs. d H3 HI [H HI). By contrast, a recent publi- 
cation [[HI concludes that the discovery of the cesium three- 
body zero-energy resonance |8j] provided the first evidence for 
the existence of Efimov's effect in nature. This assessment is 
based on a fit analysis of the measurements universally re- 
lating the resonance position of 7.2 G to the field strength of 
the minimum of ^3 just below 20 G in Fig. [3] Such a proce- 
dure relies upon the assumption that universality is preserved 
across the zero of scattering length and that all B-dependence 
of Ki(B) is included in a(B), whereas this might be strictly 
justified at most in the limit a — » +00. 



A. Three-body zero-energy resonances of 133 Cs at high field 
strengths in the vicinity of 800 G 

Besides the existing evidence for the ground state of the 
Thomas-Efimov spectrum in, for instance, 4 He molecular 
beams and ultracold gases of 133 Cs, the discovery of an excited 
state JH [Tol dH] remains an elusive goal. Measurement of the 
E2 three-body zero-energy resonance at about -9 G in Fig. [3] 
could provide an opportunity for such a discovery. Ultracold 
gases of 133 Cs prepared in the excited (F = 3, irif = -3) Zee- 
man state, however, are affecte d by inelastic atom loss due to 
both two-body spin-relaxation l27ll39ll and three-body recom- 
bination. Such a coincidence of loss mechanisms could lead 
to difficulties in the interpretation of measurements 14011 . 

For this reason, we study three-body recombination in the 
vicinity of an as yet unobserved broad diatomic zero-energy 
resonance of 133 Cs gases prepared in the (F = 3,Mf = 
+3) Zeeman ground state at comparatively high magnetic- 
field strengths on the order of 800 G. Based on the coupled- 
channels model of Ref. 12711 and a fit to Eq. ©, we obtain an 
associated singularity of the scattering length at Bo = 800.6 G, 
as illustrated in the upper panel of Fig. [5] This fit gives 
the background-scattering length and resonance width to be 
flb e = 1918 flo and AB = 83.2 G, respectively, determining 
a(B) via Eq. (fJJ. Comparatively narrow zero-energy reso- 
nances at about 735 G and 825 G in Fig. are neglected in 
this treatment. 

The lower panel of Fig. [5] illustrates the complex diatomic 
molecular energy spectrum of 133 Cs in the vicinity of 800 G. 
Dashed curves are associated with states of predominantly d- 
wave symmetry, which cause the two comparatively narrow 
zero-energy resonances in the upper panel. Solid curves refer 
to the s-wave dimer levels, which are most strongly coupled 
to the scattering continuum above the zero-energy threshold. 
Among them the 65-Feshbach molecular dimer state of en- 
ergy £b (not resolved in the lower panel of Fig. |5]l causes the 
broad diatomic zero-energy resonance whose calculated po- 
sition is 800.6 G. In order to give quantitative estimates for 
K3, we use the parametric two-channel representation of the 
.s-wave dimer levels explained in the appendix, in analogy to 
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FIG. 5: (Color online) Scattering length (upper panel) and diatomic 
bound-state energies (lower panel) of 133 Cs atoms vs. magnetic-field 
strength in the vicinity of 800 G. The solid curve in the upper panel 
refers to exact coupled-channels calculations l27ll including basis 
states of s- and rf-wave symmetries, while the dashed curve is a 
fit of Eq. l(2j to the calculated data. According to this fit, the po- 
sition of the broadest singularity of the scattering length is located 
at Bq = 800.6 G (vertical dotted line). Comparatively narrow zero- 
energy resonances near 735 G and 825 G are due to dimer levels of 
predominantly <i-wave symmetry shown in the lower panel. Their en- 
ergies with respect to the zero-energy dissociation threshold are indi- 
cated by long-dashed, dashed and dot-dashed curves. These states 
are more weakly coupled to the continuum of free atoms than 5- 
wave dimer levels. For convenience, their mutual coupling is ne- 
glected in these calculations. Solid curves in the lower panel refer to 
those 5-wave dimer levels that constitute the target states included in 
our three-body recombination calculations. The 6.s-dimer level frH\ 
causing the broadest singularity of the scattering length in the upper 
panel is unresolved on the energy scale chosen in the figure. 



our treatment of three-body recombination at low magnetic 
fields. 

Figure [6] shows calculated three-body recombination loss- 
rate constants as well as metastable Thomas-Efimov trimer 
levels causing three-body zero-energy resonances in an ultra- 
cold cesium gas at magnetic-field strengths B in the vicinity of 
800 G. Similarly to Fig. [3] these predictions refer to zero ki- 
netic energy of the three initially free colliding atoms. Based 
on such full calculations (upper panel of Fig. |6]l as well as 
a simplifying single-spin-channel treatment (lower panel of 
Fig. |6]l, we estimate the three-body zero-energy resonance as- 
sociated with the ground state (« = 1) Thomas-Efimov level 
to occur between 852 G and 859 G. Similarly to the results of 
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FIG. 6: (Color online) Calculated three-body recombination loss-rate 
constants A" 3 (upper panel) and energies of metastable trimer levels 
of 133 Cs atoms labeled E\ and E 2 (lower panel) vs. the magnetic-field 
strength B on the high-field side of the 800.6 G diatomic zero-energy 
resonance. For each magnetic-field strength zero energy in the lower 
panel coincides with the scattering threshold for three asymptotically 
free atoms. Similarly to Fig. [3] the dashed and dot-dashed curves 
labeled E d [ and E d 2 in the lower panel refer to the .s-wave dimer 
levels shown in the lower panel of Fig. [5] which constitute the target 
states included in our calculations of A* 3 . Dashed and dot-dashed 
curves in the upper panel indicate their associated contributions to 
the total three-body recombination loss-rate constant (solid curve). 
Three-body zero-energy resonance peaks are labeled by the degree of 
excitation (n = 1,2) of their associated metastable Thomas-Efimov 
trimer levels of the lower panel. 



Fig. [3] the deviations between the results of these two com- 
pletely different methods are less than 9% of the diatomic res- 
onance width AB and mainly reflect the influence of the E d l 
and E A 2 dimer levels (lower panel of Fig. [6]) on the three-body 
energy spectrum. Analogously, we estimate the position of 
the excited state (n = 2) Thomas-Efimov zero-energy reso- 
nance between 807 G and 810G. While our calculations do 
not account for the full complexity of the 133 Cs diatomic en- 
ergy spectrum, as illustrated in the lower panel of Fig. [5j we 
believe the accuracy of these estimates to be similar to those 
for the low magnetic-field strengths in Fig. [3] 

In Fig. [7] we compare the total three-body recombination 
loss-rate constants, ^3, associated with the low-field and high- 
field diatomic zero-energy resonances as functions of the scat- 
tering length, a. The overall qualitative trends of the calcu- 
lated Kj near 800 G (dashed curve) tend to be similar to those 
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FIG. 7: (Color online) Comparison between the calculated total 
three-body recombination loss-rate constants A" 3 associated with the 
low-field (solid curve) and high-field (dashed curve) diatomic zero- 
energy resonances of 133 Cs vs. the scattering length, a. Filled circles 
and squares refer to experimental data of Refs. 1 8, 28] observed at 
low magnetic-field strengths (see Fig. [3} in dilute gases of 10 nK and 
200 nK, respectively. The range of scattering lengths displayed cov- 
ers the ground-state (n = 1) and excited-state (n = 2) three-body 
zero-energy resonances of the Thomas-Efimov spectrum. The in- 
set shows an enlargement of the presently experimentally accessed 
(3, l28ll range of a(B) as well as the associated calculated and mea- 
sured A" 3 loss-rate constants. 



at low fields (solid curve). In the regime of negative scat- 
tering lengths, however, the functional dependences on a dif- 
fer considerably between the two curves, demonstrating the 
strong sensitivity of Kt, to the different ^-dependent proper- 
ties of the associated, comparatively tightly bound diatomic 
target levels. Universal behavior of ^3, which is insensitive 
to the B-dependent properties of tightly bound diatomic target 
states, occurs only in the limit of large, positive a, where both 
curves tend to agree. In this regime, three-body recombination 
is strongly dominated by the weakly bound, highest excited, 
vibrational dimer state, whose long-range wave function is de- 
termined mainly by a(B) in addition to the van der Waals dis- 
persion coefficient, Cg 11911 . These parameters also determine 
the positions of Thomas-Efimov three-body zero-energy reso- 
nances in our single-spin-channel estimates of their associated 
energy levels in the lower panels of Figs.[3]and[6][20]. 

We predict a pronounced minimum of ^3 to occur at about 
891 G in the vicinity of zero scattering length. According to 
Fig. [6] a lower estimate of the associated loss-rate constants 
is on the order of only 10~ 44 m 6 /s. Given the comparison be- 
tween theory and experiment of Fig. [3] we believe that recom- 
bination into dimer levels beyond those included in our calcu- 
lations as well as finite temperature effects could significantly 
increase the minimal value of ^3. Similarly to the experiments 
of Refs. l4ll l42ll . however, the associated minimal atom loss 
may still allow Bose-Einstein condensation of 133 Cs at high 
fields using magnetically tunable diatomic interactions. This 
possibility would be crucial to three-body recombination mea- 
surements at negative scattering lengths using the technique 



reported in Ref. H. 



B. Thomas-Efimov spectrum of 85 Rb 

Both examples of diatomic zero-energy resonances of 133 Cs 
discussed in this paper were broad and entrance-channel dom- 
inated [19]. Entrance-channel dominance allowed us to de- 
scribe and interpret near threshold three-atom resonance states 
using single-spin-channel models, largely in accordance with 
Thomas' 04J] andEfimov's [1] original treatments. Their broad 
nature should be crucial for a possible experimental resolution 
of several three-body zero-energy resonances, given present 
limitations in magnetic-field control. In both of these exam- 
ples, however, the background-scattering length was positive. 
This implies that the dominant three-body recombination peak 
refers to the previously observed ground state of the Thomas- 
Efimov spectrum, whereas the excited-state peaks occur at 
large negative scattering lengths (|a| > 20000 oq in Fig. [7} and 
accordingly large K3 loss-rate constants. Consequently, fast 
depletion of an ultracold gas at such near resonant magnetic- 
field strengths as well as thermal broadening of peaks ll43ll 
could impose limitations on the feasibility of associated mea- 
surements. Recent studies B44ll therefore suggest alternative 
measurements in mixtures of ultracold gases with unequal 
masses of the constituent particles. In such systems the num- 
ber of exotic three-body energy levels depends sensitively on 
the associated mass ratios (see Refs. II45[ 14611 and references 
therein). 
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FIG. 8: (Color online) Efimov plot associated with 85 Rb trimers I20I1 
in the vicinity of the 155 G diatomic zero-energy resonance. The 
solid and dashed curves refer to the Ej and E\ levels, respectively, 
while the dotted line indicates the energy of the Feshbach molecu- 
lar dimer state. The inset shows the magnetic-field dependence of 
the scattering length a(B) for this atomic species with a negative 
background-scattering length of a bs = -443 a l47ll . 



In accordance with Ref. [20], ultracold gases with broad, 
entrance-channel dominated diatomic zero-energy resonances 
lll9ll and negative background-scattering lengths may be best 
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suited for the observation of an excited trimer state in the 
original context of Efimov's suggestion for identical Bose 
particles of equal mass £0]. An associated Efimov plot is 
shown in Fig. [8] for the example of 85 Rb in the vicinity of 
155 G. According to the inset, the inverse scattering length 
l/a(B) increases with increasing B and extends from zero field 
(dashed vertical line) to the zero of the scattering length at 
about 165.71 G ll47l l48ll . Consequently, a range of 1/a from 
1/flbg < to +oo is experimentally accessible. As opposed 
to the cesium case in the upper panel of Fig. [4] the ground 
Thomas-Efimov level never reaches the threshold in the vicin- 
ity of Z?o = 155 G. Consequently, the only available three- 
body zero-energy resonance refers to the excited and as yet 
unobserved Ei level. 

Besides the possibility of measuring an associated three- 
body recombination loss-rate peak, the existence of the ex- 
cited Thomas-Efimov state could be verified by populating 
this trimer level using a magnetic-field sweep technique I20I1 . 
Such experiments might be performed in the gas phase as 
well as in optical lattices with three atoms occupying each 
site |2(1 SH H3, EH- Due to the larger spatial extents of 
their bound-state wave functions ll20ll . the lifetimes of excited 
£2 states of 85 Rb3, with respect to spontaneous dissociation 
via either spin relaxation lf48l l52ll or atom-dimer relaxation 
l53l I54I1 . should be significantly longer than those of ground 
states of the Thomas-Efimov spectrum. Such an increase in 
lifetime of these ultracold trimer molecules might be sufficient 
for experimental studies of their exotic properties. 



IV. CONCLUSIONS 

We have shown by an intuitive argument as well as full 
numerical calculations that the recent observation of a three- 
body zero-energy resonance in ultracold 133 Cs gases [8] refers 
to the ground level of the Thomas-Efimov spectrum that has 
been discovered in 4 He molecular-beam experiments lfl4ll . 
Unlike 4 He3, the metastable cesium trimer observed in the 
experiments of Ref. |8[] is associated with a Borromean state 
13011 . providing, to our knowledge, the first experimental ev- 
idence for their existence in molecular physics. Our studies 
show that resonance-enhanced three-body recombination in 
ultracold gases at negative as well as small positive scatter- 
ing lengths is sensitive to the strongly magnetic-field depen- 
dent Zeeman-state compositions of the comparatively tightly 
bound diatomic target states. This implies that fitting pro- 
cedures incorporating the magnetic-field dependence of the 
rate constants only through the scattering length are incom- 
plete and could therefore lead to unreliable conclusions about 
universality of trimer levels. As opposed to the magnitudes 
of the loss-rate constants, the mere magnetic-field positions 
of Thomas-Efimov zero-energy resonances largely follow the 
qualitative trends implied by three-body universality, irrespec- 
tive of the details of diatomic energy levels far away from the 
scattering threshold. 

The long-standing problem of detecting excited Thomas- 
Efimov states ifTol [ill [l5tl remains an elusive goal. Based 
on numerical predictions of three-body recombination loss- 



rate constants, we have investigated the possibility of observ- 
ing associated three-body zero-energy resonances in ultracold 
133 Cs gases. Our studies focused on a broad, entrance-channel 
dominated [ 19] diatomic zero-energy resonance, which is pre- 
dicted to occur near 800 G. In this case, the background- 
scattering length is positive. We have shown that for this 
reason an associated experiment would require observation of 
two distinct three-body recombination loss-rate peaks. The 
excited state three-body zero-energy resonance occurs at large 
negative scattering lengths. Thermal broadening [43] as well 
as strong atom loss could therefore impose limitations on 
the feasibility of associated measurements. We reach the 
conclusion that ultracold gases involving entrance-channel 
dominated diatomic zero-energy resonances with negative 
background-scattering lengths, such as 85 Rb, may be best 
suited for the first observation of excited Thomas-Efimov 
trimers. Their detection might be achieved by populating such 
metastable states, with sufficient lifetimes, via magnetic-field 
sweeps [20] in ultracold gases or in optical lattices with three 
atoms occupying each site I20I l49l [501 l5lll . 
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APPENDIX A: AGS APPROACH TO THREE-BODY 
RECOMBINATION 

1. Scattering matrix 

a. Asymptotic scattering channels 

Throughout this appendix, we consider collisions between 
three identical Bose atoms. Three-body scattering is charac- 
terized by the arrangements of the atoms on asymptotically 
large time scales before and after a collision. These arrange- 
ments determine the asymptotic entrance- and exit-scattering 
channels. Such channels can involve a bound molecular pair 
of atoms asymptotically spatially separated from the third 
atom or three free atoms. In the case of three-body recombina- 
tion the entrance channel consists of three asymptotically free 
atoms, whereas the exit channel contains a diatomic molecule 
and a free atom. 

To characterize the positions of the atoms in the different 
asymptotic scattering channels, we choose the Jacobi coordi- 
nates of Fig. [9] in the center-of-mass frame. These coordi- 
nates are particularly well suited to describe asymptotic ar- 
rangements involving a bound pair and a free atom. In the 
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FIG. 9: Jacobi coordinates of the relative motion of three atoms. The 
set of coordinates p = p { and r = is selected in such a way that 
it describes the hypothetical situation of an interacting pair of atoms 
(2, 3), whereas atom 1 is considered as a spectator. Two other sets of 
these coordinates describing different arrangements of the atoms can 
be obtained through cyclic permutations of the atomic indices. 



following, we label the sets of coordinates by the Greek index 
a = 1,2,3 of the free atom. For convenience, the particular 
set of coordinates shown in Fig. [9] will be denoted, without 
explicit reference to its index, simply by (p, r). 



b. Transition matrices 

We assume the three atoms to interact pairwise depending 
on their internal Zeeman states [23]. Consequently, the three- 
body Hamiltonian in the barycentric frame, 



H = H + Vi + V 2 + V 3 , 



is comprised of a kinetic energy contribution, 



H = |bg) 



h 2 , h 2 
v 2 

<2 m \ P 



2(|m) 



2(f) 



<bc 



(Al) 



(A2) 



in addition to the pair potentials Vu Vi, and V3. Here m is 
the atomic mass, |bg) denotes the entrance-channel spin prod- 
uct state, consisting of the three one-body Zeeman states in 
which the atoms of the gas are prepared, and V\ is the poten- 
tial of the atom pair (2, 3). The remaining interactions V2 and 
V3 are determined through cyclic permutations of the atomic 
indices. Each potential couples entrance- and closed-channel 
Zeeman state configurations in Feshbach resonance-enhanced 
collisions B23I1 . As the entire physics described in this pa- 
per involves energy ranges in the close vicinity of the scat- 
tering threshold, and accordingly large interparticle distances, 
we have neglected short-ranged, genuinely three-body forces 
in our approach ifTTi EBll . 

The probabilities for transitions between the different 
asymptotic scattering channels can be conveniently described 



in terms of the scattering matrix S. To determine the S- 
matrix, we follow the approach by Alt, Grassberger and Sand- 
has (AGS) f56ll57ll and introduce the following Green's func- 
tions associated with the different asymptotic arrangements of 
the atoms: 



G (z) = (z-H r\ 
G a {z) = (z - Ho - V a y 



(A3) 
(A4) 



Here z is a complex variable with the dimension of an en- 
ergy. Throughout this appendix this variable is given by a 
limit, z = Ei + iO, indicating that the energy of the three ini- 
tially noninteracting atoms, Ei, is approached from the upper 
half of the complex plane. The free Green's function Go(z) is 
associated with three asymptotically free atoms, while G a (z) 
describes the configuration in which the atom with the index a 
is free and the remaining atom pair interacts via the potential 
V a . The complete Green's function 



G(z) = (z-fl)" 



(A5) 



implicitly determines the AGS transition matrix U a p(z) 
through the relation: 



G(z) = 6 afs Gp(z) + G a (z)U a p(z)Gi3{z). 



(A6) 



Here the indices (3 — and a — are associated with entrance 
and exit channels involving three free atoms, respectively. The 
indices /?, a = 1,2,3 correspond to arrangements in which 
the atom with the index /3 or a is free, respectively, while the 
remaining atom pair is bound to a molecule. 

The transition matrices U a p(z) determine the elements of 
the S -matrix associated with the different asymptotic scatter- 
ing channels. The general derivations of Refs. |p6[|57fl show, 
in particular, that the S -matrix element for three-body recom- 
bination into a diatomic molecular bound state |0d) is given by 
the formula: 



ff <<jf, (f> d \S Iqi, pi, bg> = - 2m5(Ef - Ei) 

x t ,<qf,<Mt4o(z)lqi,Pi,bg>. 



(A7) 



Here |q;,pi,bg) = |qi)|pi)|bg) and |q f , d >a = lqf>al^d>« 
are the general unsymmetrized initial and final momentum- 
normalized energy states associated with the entrance and exit 
channels, respectively. Consequently, p and q denote the mo- 
menta associated with the Jacobi coordinates r and p, respec- 
tively. This implies, for instance, that qf is the momentum 
of the center of mass of the molecule produced relative to the 
free atom with index a, and (p\qt) = exp(/qf • p/S)/(27r^) 3 ^ 2 
is the associated plane wave. Similarly, \(f>d) a is the diatomic 
molecular state of the interacting atom pair in the exit chan- 
nel. We have arbitrarily chosen the Jacobi momenta describ- 
ing the entrance-channel state in accordance with the selection 
of atom 1 in Fig.|9]and omitted the index associated with this 
particular arrangement of particles. The asymptotic energies 
of the three atoms in the barycentric frame are symmetric with 
respect to permutations of particles and related to the Jacobi 
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momenta by the formulae: 



Ei 



Am 

E - 3 A 

4-m 



Pi 



(A8) 
(A9) 



Here E& is the binding energy associated with the di- 
atomic molecular target state \<f>&) a . Consequently, the 
entrance- and exit-channel states solve the Schrodinger equa- 
tions ff |qi, pi.bg) = £i|q;, pi.bg) and (H + V a )\q f ,(f>d) a = 
Ef |qf, <pd)a> respectively. The S -matrix elements for transi- 
tions between other scattering channels can be represented in 
terms of the transition matrices of Eq. (IA6t through formulae 
similar to Eq. dA7} 01 HI. 



2. Three-body recombination loss-rate constant 

a. Probability for the recombination of three atoms 

To determine the three-body recombination loss-rate con- 
stant, we consider three initially free atoms in a large box of 
volume *V which will later be taken to be infinite. The prob- 
ability Pfi for recombination into the molecular bound state 
\<pd), in a time interval of length Af, can be obtained using 
Fermi's golden rule. This yields: 



P& = 



18 (2nh) 6 
¥ <V 2 



I 



dq f 



'sin[i(£ f -£i)A?/fcp 



\{E f - Em 

x |i<q f ,0 d |t/i,o(z)S|qi,Pi,bg)| 2 . 



(A10) 



We note that the complex argument "z = E\ + iO" of the 
transition matrix E/^oCz) indicates that the physical energy 
Ei is approached from the upper half of the complex plane. 
To account for the identical nature of the Bose atoms, the 
entrance-channel state |qi,p;,bg)i(2^) 3 / r V in Eq. ( IAIOI i is 
symmetrized by the projection 



Is* 



(All) 



where the sum extends over all 3 ! =6 permutations f of the 
atoms. Due to the symmetric nature of the projected initial 
wave function in Eq. flAlOt the index a — 1 associated with 
the final atomic arrangement is arbitrary. 

The approximations underlying the Fermi golden rule argu- 
ment and the associated continuum limit of the energy states 
of the box determine the range of validity of Eq. flAlOl i. The 
time scale At is limited from above by the requirement that 
phenomena associated with secondary collisions of atoms due 
to reflection from the boundaries of the box are negligible. 
This implies that At must be much smaller than the ground- 
state period of the box. This gives the restriction 



where E^° x = tfn 2 /(2m r V 2 ^ 3 ) is the single-particle ground- 
state energy of the box. Given that ultracold collision energies 
Ei are usually far below the breakup threshold, i.e. Ei <k \E&\, 
a simple calculation using Eqs. (IA8b and jA9i shows that the 
condition 



%l |£ d | «: At 
allows us to perform the formal limit 

'sin^iEf- Ei)At/tiW~ 



(A13) 



(Ef - Eo/n 



2nftd(Ef - Ei)At (A14) 



in the integral over the final momentum qf of Eq. ( lAlOt . The 
applicability of Eq. (IA14I) ensures that the transition proba- 
bility per unit time, i.e. Pn/At, is independent of Af within 
the limitations set by Eqs. dA12b and dA13b . We note that the 
lower bound on Af set by Eq. ( IA1 3b is also necessary for the 
validity of Eq. dAlOl ). The time scale %j \E&\ is therefore char- 
acteristic for a three-body recombination event. 



b. Three-body recombination loss rate for a thermal Bose gas 

The interpretation of atom loss observed in a dilute atomic 
gas in terms of rate phenomena associated with any chemi- 
cal reaction crucially relies upon the possibility of introduc- 
ing coarse-grained scales in space and time. To this end, we 
divide the volume of the gas into regions of virtually constant 
density n(x) characterized by their central position x. The vol- 
ume *V of each region sets the scale for the unit of length in 
the coarse-graining approach. The unit of time Af is limited 
by Eqs. (IA12b and ( IA1 3b except that the volume of the box in 
Eq. dA12t needs to be replaced by [n(x)] . This replacement 
rules out collisions of three specific colliding atoms with any 
other constituents of the gas during a time interval of length 
Af. Equations (IA12b and (IA1 3b thus determine the condition 
for the applicability of a rate treatment of three-body recom- 
bination in a dilute gas to be: 

\E d \ » [n(x)] 2/3 H 2 /tn. (A15) 

In the absence of diatomic zero-energy resonances Eq. (IA15b 
is fulfilled for most realistic sets of physical parameters of di- 
lute ultracold gases. Magnetic tuning of the bound-state en- 
ergy £b of the highest excited vibrational state, however, may 
violate this condition. 

Given that Eq. (IA15b applies, we follow the treatment of 
Ref . 115811 to derive the loss-rate constant for three-body recom- 
bination. This treatment relies upon a summation of the mi- 
croscopic transition probabilities of Eq. (lAlOb . per unit time 
Af, over the number a N 3 /6 of atomic triplets in a re- 
gion with atoms. This yields the number of recombination 
events. Under the assumption that all three atoms are lost in 
each recombination event, the associated loss-rate constant in 
a thermal ultracold Bose gas is given by the formula: 



^3 



i^Pfi/Af 
12^m 



Af « h/E b o 0X = 2nfV 2li l(n 2 h), 



(A12) 



(2nh) 6 q t J rfQ f |i<q f ,0 d |E/i, o (z)|O>| 2 . (A16) 
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Here dClt denotes the angular component of d(\f . In the deriva- 
tion of Eq. dA16t we have taken advantage of the fact that 
the transition amplitude i(qf, 0dl^i,o(-Ei + /0)5|q;,pi,bg) of 
Eq. dAlOt is often insensitive with respect to typical collision 
energies E\ of a cold gas. For convenience, we have thus ap- 
proximated the entrance-channel state by the zero-momentum 
plane wave |0) of the relative motion of three noninteracting 
entrance-channel atoms and performed the limit z = ;0. The 
energy conservation implied by Eq. ( IA14I ) then determines the 
modulus of the final momentum of the diatomic molecule rel- 
ative to the third, free, atom to be: 

qt = 2 Vm|£dl/3. (A17) 

The rate equation for the total number of atoms in a cold gas, 
N(t), can be obtained through averaging the local rate equa- 
tions over all regions of constant density. This local density 
approach yields Eq. (|T). If recombination events occur in 
a Bose-Einstein condensate then the three-body recombina- 
tion rate K^{n 2 (t)) needs to be divided by a factor of 3! = 6 
ll59il60ll . because the zero-energy entrance-channel state | 0) i s 
already symmetric with respect to permutations of atoms 115811 . 



3. AGS method 

a. AGS equations for transition matrices 

In the following, we describe the general approach we have 
employed to determine the transition amplitude of Eq. ( |A16t 
based on the AGS method lf56T,l57ll . Other techniques for solv- 
ing the three-body Schrodinger equation in the context of ul- 
tracold collisions often rely upon h ypersph erical-coordinate 
representations (see, e.g., Refs. l6lp7l l6lll62ll63Tl ). According 
to Refs. 11561 15711 . all complete transition matrices fulfill the 
coupled AGS equations: 

3 

U ap (z) = (1 - 6 a p)G \z) + J] U ay (z)G y (z)V y , (A18) 

7=1 

rW 
3 

U ap {z) = (1 - 5 ap )G Q l {z) + J] V y G y {z)U yP {z). (A19) 

7=1 

Here both Eq. ( IA18b and Eq. ( IA19b determine the transition 
matrices U a p(z) with indices a,/3 = 1,2,3. We shall refer 
to them as the prior and post versions of the AGS equations, 
respectively. Given the solution to these equations we insert 
P = into the prior version of Eq. dA18l ). This yields: 

3 

Uao(z) = G \z) + J] U ay (z)G y (z)V y . (A20) 
y=i 

Taking advantage of the invariance of the zero-energy plane- 
wave state with respect to permutations of the atoms, the AGS 
equation ( IA20b gives the transition amplitude of Eq. ( IA16b to 



be: 

3 

i<q f , 0d|C/i,o(z)|O> = Yj i<«Jf. UU,p(z)G p (z)Vp\0) 

13=1 
3 

= J] „<q f , </> d \U al (z)G l (z)V 1 \0). (A21) 

a=l 

We note that each transition matrix U a p(z) is equivalent to 
the complete three-body Green's function due to Eq. ( IA6b . 
The solutions of the coupled sets of equations ( IA1 8b as well 
as jA19b thus determine all energy levels of the three-body 
Hamiltonian rather than just the transition amplitudes between 
asymptotic channel states. 

b. AGS equations for the three-body recombination transition 
amplitude 

To render the problem of determining the exact amplitude 
of Eq. ( IA2 1 b into a more practical form, we introduce the fol- 
lowing wave functions: 

\*ai) = t4i(z)Gi(z)Vi|0>, (A22) 
l<ALi> = G Q \z)G l {z)V l \Q){\ - S al ). (A23) 

Here the index a assumes the values a - 1,2,3. A straight- 
forward calculation using Eq. jA19b shows that \\jj a \) is deter- 
mined by the coupled set of Faddeev-type II 6411 equations: 

3 

Ifcrl) = > + Yj VyG 7 (zMyl)- (A24) 

7=1 
7*0 

The solution of Eq. ( IA24I ) gives the transition amplitude for 
three-body recombination to be: 

3 

i<q f ,0d|f/i,o(z)|O> = 2«<«lf'^l^i>. (A25) 

(7=1 

We use Eq. (IA24b to determine the associated loss-rate con- 
stant of Eq. dA16b . the implementation of which we discuss in 
Subsection [A3] In the following subsection, we first provide 
a detailed description of the interatomic potentials. 

4. Resonance-enhanced three-body recombination 

a. Three-body spin channels 

Ultracold resonance-enhanced diatomic interactions lfl6l 
El d usually occur due to coupling between the entrance- 
spin-channel and a single energetically closed channel. For all 
the broad diatomic zero-energy resonances described in this 
paper the physical origin of this coupling is predominantly 
spin exchange lfT7il65ll . In the context of three -body collisions 
the resonant closed channel is characterized by a Zeeman- 
state configuration, |cl)„. Here the index a - 1,2,3 indicates 
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that atom a is in the Zeeman state in which the Bose gas is pre- 
pared, whereas the remnant atom pair is in the closed-channel 
spin state, whose strong coupling to the entrance channel 
causes a diatomic zero-energy resonance 

Emm. The 

three-body spin state is orthogonal to the product state |bg), 
i.e. tl .(cl|bg) = 0. In the case of, e.g., the 133 Cs low-field 
zero-energy resonance of Figs. [2]and[3]the resonant diatomic 
closed channel is comprised of Zeeman states mainly from the 
upper F = 4 hyperfine level, whereas the atoms are prepared 
in the (F = 3, m F = 3) Zeeman ground state. This implies the 
orthogonality relation, 



a {cl\c\)p = 5 ap , 



(A26) 



for all indices a,B — 1, 2, 3 of the three-body closed-channel 
spin states. We assume in the following that the resonant inter- 
channel coupling fulfills Eq. dA26l ), i.e. no single-atom Zee- 
man state is shared between the diatomic entrance and closed 
channels. This assumption does not affect our description of 
diatomic energy levels for which the closed channel is always 
orthogonal to the entrance-spin-channel, even when single- 
particle Zeeman states are shared between them. It facilitates, 
however, our calculations of three-body recombination loss- 
rate constants in 133 Cs gases and is strictly valid only for the 
low-field resonance shown in Figs.|2]and[3] 



b. Single-resonance approach 

Resonant enhancement of ultracold collisions is often well 
described in terms of near degeneracy of a single diatomic 
closed-channel energy state, Ifoes), with the dissociation 
threshold of the entrance channel lfm[T7l,[l9tl . In an extension 
of this concept to three-body scattering the resonance state, 
l^res, cl) a , refers to the situation in which a pair of atoms occu- 
pies the state \<f> ies ), whereas the third atom with index a plays 
the role of a spectator. Throughout this appendix we choose 
|^ re s) to be unit normalized, i.e. (0 re sl0res) = 1. In accordance 
with such a single-resonance approach, the pair potentials of 
the three-body Hamiltonian ( IAU are of the following general 
formlil: 

V a =|bg>^ g <bg| + |0res,cl> a £resCB) 
res? v ^ 1 l 

+ lAc, cl> a „<0 res , bg\W a . (A27) 

Here V^ 8 is the background-scattering potential, W a (r a ) is the 
potential energy of interchannel coupling depending on the 
relative position r a of atoms of an interacting pair, and E res {B) 
is the resonance-state energy [19]. As we have chosen zero en- 
ergy to coincide with the scattering threshold for three asymp- 
totically free atoms, the three-body Hamiltonian ( IA1 b depends 
on the magnetic-field strength only through E K!i (B). To a good 
approximation this dependence is determined by 12511 : 



Eres(B) - /i r es(fi ~ B les ). 



(A28) 



Here yu les is the difference in magnetic moment of the reso- 
nance state and a pair of asymptotically separated entrance- 
channel atoms, and B les is the crossing point between E res and 
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FIG. 10: (Color online) Comparison between coupled-channels 
(filled circles) and two-channel single-resonance (long-dashed and 
dot-dashed curves) approaches to the highest excited vibrational s- 
wave states of 133 Cs 2 vs. the magnetic-field strength B in the vicinity 
of 800 G. The dashed line labeled £ les (B) indicates the bare Feshbach 
resonance-state energy, while the horizontal dashed line refers to 
the second highest excited bare vibrational level of the background- 
scattering potential, £ b |. The energies of the associated highest ex- 
cited bare level, E bg v and of the Feshbach molecular dimer state, 
£b(S), are not resolved. The point of degeneracy of £b with the 
scattering threshold (zero energy) is indicated by the vertical dotted 
line labeled So. It coincides with the position of the broad singular- 
ity of the scattering length in Fig. [5] This measurable zero-energy 
resonance position is distinct from the zero-energy crossing point 
of the bare Feshbach-resonance level, S res . The deviation at lower 
B-fields between the coupled-channels and two-channel E d 2 levels is 
mainly due to a strong avoided crossing of dressed, coupled-channels 
.s-states in the vicinity of 600 G, which is not included in the two- 
channel approach. 



the entrance-channel scattering threshold, i.e. £ , res (B les ) = 0. 
Figure [10] shows that for the entrance-channel dominated di- 
atomic zero-energy resonances [ 19] studied in this paper, B res 
is considerably shifted with respect to the measurable position 
Bo of the singularity of a(B) referred to in Eq. ©. 



5. Implementation of the AGS equations 



a. Interchannel coupling 



For convenience, the contributions to Eq. ( lA27b describing 
interchannel coupling can always be represented in terms of 
an amplitude and a wave function \x), i.e. 



w|0 res > = \xX- 



(A29) 



Within a limited range of energies close to the dissociation 
threshold the precise functional form of x(f) - (r\x), asso- 
ciated with its dependence on the distance r = |r| between 
atoms of an interacting pair, is unresolved. As a consequence, 
several different models of coupling are capable of predicting 
the same low-energy physics. For this reason, it is sufficient 
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for all applications in this paper to describe the interchannel 
coupling simply in terms of an amplitude and a range. In ac- 
cordance with Ref. i66ll . we choose the following arbitrary but 
convenient Gaussian function in momentum space: 

- P V/(2S J ) 

XV = <Pl*> = ^JT- (A30) 

Here <x is the associated range parameter, p denotes the mo- 
mentum of the relative motion of an atom pair, and p = |p| is 
its modulus. The model parameters £ and <x are determined 
mainly by the resonance width AB of Eq. (f2j) and the reso- 
nance shift Bo - fi res |fl9jl. Interchannel coupling is also re- 
sponsible for the avoided crossing of the levels E d { and E_ 2 , 
illustrated in Figs. |2l and [TOl 

b. Background scattering 

The background-scattering potential Vbg describes the in- 
teraction of a pair of entrance-channel atoms in the absence 
of interchannel coupling at magnetic fields asymptotically far 
away from diatomic zero-energy resonances. In combination 
with the resonance state, its bare vibrational energy levels 
determine the number and main properties of comparatively 
tightly bound dressed, two-channel target states for three- 
body recombination at negative a. Their significance tends 
to decrease with increasing binding energy. Whereas realistic 
alkali atom pair potentials typically support dozens of dimer 
levels, we found it sufficient to account for only the two high- 
est excited vibrational states in our model of Vb g to determine 
the correct magnitudes of K^(B) in Fig. [3] To this end it is cru- 
cial, however, to precisely mimic the properties of the dressed 
dimer energy levels £ d 1 and E d 2 indicated in Figs. |2l and [TOl 
which have been inferred from exact coupled-channels calcu- 
lations ||27ll . 

Similarly to the interchannel coupling, the precise func- 
tional form of the microscopic potential Vb g (r) is unresolved 
for the energy range of interest, leaving a choice of models 
predicting the same physics. To describe background scat- 
tering as well as the two highest excited vibrational levels 
of Vbe(r) we employ the following separable representation 
ll67l 16811 of the background-scattering potential: 

2 

V* = J>y>£/<*yl, (A31) 
7=1 

Such pseudo-interactions are particularly useful in numeri- 
cally exact treatments of the three-body Schrodinger equation 
and can be strictly motivated by complete separable expan- 
sions of diatomic transition matrices l57l l69ll . In accordance 
with Ref. B66I1 . we choose each form factor \xj) of Eq. ( 1A3U 
to be a Gaussian function in momentum space: 

e - P V;/(2s 2 ) 

xAp) = <Pl*;> = (2nn)V2 ■ (A32) 

The amplitudes and range parameters crj are determined for 
j = 1 , 2 in such a way that Eq. (IA3U recovers a\, g of Eq. (O, 



the effective range of the microscopic background-scattering 
potential I70t 17111 . as well as its two highest excited vibrational 
energy levels. 

Given the resonance state energy E res (B) of Eq. dA28l >. 
Eqs. dA29), dA30l . dA3Tb and dA32l completely set up our 
single-resonance model of the pair interaction ( IA27b . We note 
that the associated diatomic two-channel bound-state energies 
do not depend on the functional form of (f> K s(r) II 1911 . whereas, 
in general, calculations of K3 loss-rate constants could be sen- 
sitive to it. It turns out that the orthogonality assumption of 
Eq. ( |A261 > allows us to treat three-body recombination with- 
out specifying the wave function ^resM- Similar phenom- 
ena have been discussed in the context of pairing in ultracold 
Fermi gases IT72ll . The predictive power of our model has so 
far been demonstrated not only in applications to two-body 
physics 01911 . but also to resonance-enhanced three-body de- 
cay of 87 Rb Bose-Einstein condensates 1I23I1 as well as atom- 
dimer collisions In contrast to the broad diatomic 133 Cs 
and 85 Rb resonances of the present paper, the narrow 87 Rb 
1007G resonance of Refs. 12311 and 17311 is closed-channel 
dominated fl9Tl . 

c. AGS integral equations 

Given the pair interactions of the three-body Hamiltonian 
(IAU . the general AGS equations (IA24b can be converted into 
integral equations by multiplying them by a (q, p| from the left 
and inserting the completeness relation, 

J dp Jrfp|q,p> 77 <q,p| = 1, (A33) 

for each y + a into the sum on the right-hand side. Using a 
partial wave expansion, the numerical solution of Eq. (IA241 > 
is equivalent to matrix inversion on a two-dimensional grid 
involving both Jacobi momenta q = |q| and p = |p|. This pro- 
cedure is commonly known as the momentum space Faddeev 
approach B7411 . 

Similarly to the usual separable-potential approach f57l 
l69ll . it is the particular form of Eq. dA27t in combination 
with Eq. (IA26b that allow us to reduce this problem to a one- 
dimensional integral equation in the variable q. To this end, 
we introduce the following complete transition amplitudes, 

Xbg,;(q) = ~2^ a {q,Xj>bg\G a (z)\i// aA ), 

1^ 

X d (q) = - 2j a{q, ^ res , cl|G„(z)|^i>, (A34) 

a=l 

and inhomogeneous terms, 

*ci(q) = ~J]M,<r>™M\G a (zW a<l ). (A35) 

a=l 
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Here the summation over the index a = 1,2,3 refers to the 
equivalence of particle arrangements due to the identical na- 
ture of the three Bose atoms. Given these amplitudes, mul- 
tiplying the AGS equations dA24t by a {q,Xj,bg\G a (z) and 
a (q,<p KS ,cl\G a (z) from the left and summing over the index 
a yields the following integral equations: 

2 r 

k=i J 

+ j dq"K^ c \q,q')X cl (q'l 

X c i(q)=l 1 cl (q) + 2 f dq"K?' bg (q,q')X bg , k (q') 
k=i J 

+ J" rfq'#ci,ci(q,q')*ci(q'). (A36) 

Here the index j assumes the values 1 and 2 and the integral 
kernels are given by the following expressions: 

*?<J' bg (q,q') =2 2 <q,^,bg|G 2 (z)|q',^,bg) 1 ft, 
*7? g - d (q,q') =2 2 <q,^,bg|G 2 (z)|q',^,bg) 1 ft 
^(q, q') =2 2 <q, res , cl|G 2 (z)|q',^, bgh ft, 
*iw(q, q') =2 2 <q, res , cl|G 2 (z)|q',#, bg), ft (A37) 

The indices a = 2 and y = 1 of the final and initial particle 
arrangements in these integral kernels are arbitrary as long as 
they differ, and the prefactors of 2 stem from the summation 
over indices y + a in the general AGS equations dA24l >. As 
the kernels ( IA37I ) do not couple different partial waves associ- 
ated with q, Eqs. ( IA36t can be treated, in general, as indepen- 
dent sets of integral equations in the variable q = \q\ for each 
partial wave individually. Since we have restricted ourselves 
to zero-energy initial conditions (z = iO) the inhomogeneous 
terms are isotropic, i.e. ^ g /q) = ^ g ,/?) and Xjj(q) = X^ (q). 
Consequently, for all applications in this paper Eqs. JA36I ) in- 
volve the i-wave only. 



d. Numerical treatment of the AGS integral equations 

Due to the arrangement-channel Green's function, G 2 (z), 
the kernels of Eqs. (IA37t have poles at the physical ener- 
gies associated with dressed, coupled-channels bound states, 
as indicated in the upper panel of Fig. [2] In order to explic- 
itly locate the positions of these singularities, we introduce 
the effective diatomic energy variable of the relative motion, 
E = z — 3q 2 /(4m), as well as the entrance-channel Green's 
function, 



C bg (E) = {E + H 2 V 2 /m - V hg y 



(A38) 



According to Ref. B66H . the resonance denominator of the 
complete diatomic Green's function G 2 (z) is given by 



At each magnetic-field strength B, its zeros determine the 
locations of all diatomic bound-state energies E& = E^{B) 
through S(B, Ed) = 0. The dressed, coupled-channels dimer 
states are characterized by the wave-function normalization 
coefficient, 



Z(B, E) = 



1 



(A40) 



,1E 



which determines their closed-channel admixtures in the 
lower panel of Fig. [2] via the relation Z(B) = Z(B, E&). 
This coefficient in combination with the arrangement-channel 
Green's function and its resonance denominator determine the 
final state of Eq. (IA251 > to be 

lq, <Pi) a = G a {z)\q, res , cl>„ S(B, E) yjZ{B,E) (A41) 

in the limit E — > Ed(B)- 

Besides the poles in Eqs. ( IA35b and ( IA371 > at the diatomic 
bound-state energies, the resonance denominator of Eq. JA39b 
as well as the coefficient of Eq. (lA40b involve singularities 
due to the bare Green's function of Eq. JA38I ). Although these 
unphysical singularities associated with bare vibrational levels 
cancel exactly in the AGS integral equations ( IA36b . for any 
numerical implementation of Eqs. JA36b it is imperative to 
explicitly treat these cancellations. To this end, we introduce 
the modified, dimensionless transition amplitudes, 

X bgJ (q) = qX bgJ (q)^^= ^Z{B,E), 
X d (q) = qX d (q)S(B, E) yjZ(B,E), (A42) 
in addition to the following inhomogeneous terms: 

KM) = qKMW- E ) jz(B,E). (A43) 

Here the index j assumes the values 1 and 2, similarly to 
Eqs. ( IA34I ) and dA351 >. The associated AGS integral equations 
can be readily obtained from Eqs. JA36I ), which yields: 

X bgJ (q) =Xl g j (q) + 2/ dq' K)f° % {q, q')X HJc (q') 

k=l ^° 

+ dq"K^' c \q,q')X cl (q'), 
Jo 

cl (q) =X i cl (q) + V dq' flf bg ( 9 , q')X bgM {q') 

+ dq"K a , ci (q,q')X ci (q'). (A44) 
Jo 



X, 



6(B, E) = E- E KS (B) - |<r|%|G bg (£)^>. 



(A39) Here the modified kernels are given in terms of Eqs. ( IA37I ) by 
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the following expressions: 



bg , b , Z(B,E) EtB,E) <Xk\Xk 



<Kf c \q,q') =qq'. 



Z(B,E) 8(B,E) £ 
yjz(B,E')E(B,E') V<^> 



j dn"Kf' e \q,q'), 



1>b , Z(B,E) 8(B,E) r—— 

% (g>g)=gg^ z(B>£0fi(B>£0 ftV(y^> 



'^, d ' bg (q,q'), 



VZ(fl,£') 8(B,E')^ 
J dCl"k cld {q,q')- 



(A45) 



Here c/Q' denotes the angular component of dq' and stems 
from the momentum integrals of Eqs. ( IA36t . 

The transformed one-dimensional AGS integral equa- 
tions dA44l ) are numerically stable and can be readily solved 
for the complete transition amplitudes of Eqs. (1A421 >. Using 
Eqs. dA16t and dA25l ), their solution determines the three- 
body recombination loss-rate constant via the following re- 
lation: 



K 3 = \%Tt\2nKfm% 



X cl (q t ) 



(A46) 



Here the final Jacobi momentum qi is given by Eq. ( IA17I ). 
Our implementation of the AGS method was used to deter- 
mine the magnetic-field dependence of the Kt, loss-rate con- 
stants of Ref. B23I1 as well as in the upper panels of Figs.[3]and 
[6] A similar treatment using the same kernels but different 
inhomogeneous terms has been used to predict the scattering 
amplitudes for atom-dimer collisions in ultracold Bose gases 
of 87 Rbl|2l. 



6. Parameters for 133 Cs Feshbach resonances 



a. Diatomic entrance and closed channels 



spatially homogeneous nature of the magnetic field, all pair 
potentials are invariant with respect to rotations about the field 
axis. Both the high- and low-field 133 Cs Feshbach resonances 
studied in this paper interact with the diatomic entrance-spin- 
channel mainly via spin exchange lfl7l l65ll . These interac- 
tions do not couple different partial waves associated with the 
relative motion of an atom pair. Consequently, the entrance- 
spin-channel, {F\,m\\F2,mi) = (3, 3; 3, 3), associated with 
a pair of 133 Cs atoms 1 and 2 is coupled only to those Zee- 
man states which conserve the sum of single-particle angular 
momentum projection quantum numbers, m\ + m.2 — 6. All 
possible strongly coupled closed channels are therefore char- 
acterized by the pairs of atomic quantum numbers (3, 2; 4, 4), 
(3, 3; 4, 3), (4, 3; 4, 3), and (4, 4; 4, 2). Their associated mag- 
netic moments can be estimated by the Breit-Rabi formula 
IT76ll to determine the bare resonance-state energies E Tes (B) of 



Eq. (|A28]l [19]. 



b. Low-field zero-energy resonance 

The low-field zero-energy resonance at the negative 
magnetic-field strength j26ll of Bo = -11.2G is caused by 
a Feshbach resonance level, E les (B), comprised of Zeeman 
states mainly from the excited F = 4 hyperfine level. In- 
dependently of the specific spin composition of the associ- 
ated strongly coupled closed-spin-channel state the Breit-Rabi 
formula yields yu res = h x 4.2 MHz/G for the magnetic mo- 
ment difference of Eq. (IA28t . Coupled-channels scattering 
calculations of the resonance-enhanced scattering length a(B) 
including basis states of both s-wave and af-wave symme- 
tries determine the parameters a bg and AB given in Subsec- 
tion [HA] Besides a be , the energies of the two highest excited 
bare vibrational entrance-channel states of —h xO.Ol MHz and 
—h X 1 1 1 MHz, as well as the effective range r e ff = 252 ao 
determine the separable background-scattering potential of 
Eq. ( IA3U . Given the background-scatterin g le ngth, both 
of the bare energies [70] as well as r e ff |p70t l7lll are deter- 
mined mainly by the van der Waals dispersion coefficient 
C(, = 6890 a.u. of cesium [39] (the atomic unit of Ce is 1 a.u.= 
9.5734 x 10~ 26 J nnA The resonance width AB and the shift 
Bq - B les - -29 G III9I1 determine our model of the interchan- 
nel coupling. All these physical parameters give the model 
parameters of Eqs. dA29l , dA30l >. ( TATTi and dA32l to be: 
£1 = -5.4187 x4n V2 h 2 cri/m, & = -12.9603 x 4-7t 3/2 h 2 o- 2 /m, 
£ = 0.59623807 x 2n V4 n 2 /(mcr 1/2 ), <r x = 82.3601 a , cr 2 = 
36.4963 flo, and cr = 14.4694 a . 



The nuclear spin quantum number / = 7/2 of I33 Cs IT75ll 
gives rise to two hyperfine manifolds labeled by the total 
atomic angular momentum quantum numbers, F — 3 and 
F = 4. Due to the magnetic field, each hyperfine level is 
split into Zeeman sublevels. Interatomic collisions in an ul- 
tracold gas are sensitive to the Zeeman state in which the 
atoms are prepared. Throughout this paper, this Zeeman state 
is characterized by the (F = 3,nif - +3) angular momen- 
tum quantum numbers of the atomic state which it adiabat- 
ically correlates with in the limit of zero field. Due to the 



c. High-field zero-energy resonance 

The high-field zero-energy resonance at Bo = 800 G ex- 
ists due to strong coupling between the diatomic (3, 3; 3, 3) 
entrance-spin-channel and a closed spin channel comprised 
of the (3, 2; 4, 4) and (3, 3; 4, 3) Zeeman-state configurations. 
According to coupled-channels calculations of the s-wave 
bound states, their relative admixtures to the closed-channel 
superposition state are about 57% and 43%, respectively. 
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Based on the Breit-Rabi formula, these admixtures imply an 
estimated magnetic-moment difference associated with the 
Feshbach-resonance state of yU les = h x 2.6 MHz/G. We note 
that due to the single-particle spin overlap between the en- 
trance and closed channels, Eq. ( lA26b does not strictly ap- 
ply to the 800 G high-field zero-energy resonance. We be- 
lieve, however, that the associated sensitivity to short-distance 
physics, such as the functional form of the resonance state, 
<pies(r), is sufficiently small for our calculations to predict the 
correct orders of magnitude of ^3. Our estimates of three- 
body zero-energy resonance positions are unaffected by spin- 
overlap, as they only involve energy scales much closer to 
the dissociation threshold than three-body recombination into 
tightly bound diatomic target states at negative a(B). 



The physical parameters a\, s and AB determining the 
resonance-enhanced scattering length of Eq. © near 800 G 
are given in Subsection IIII A I Based on a\, g and C& we 
calculate the energies of the two highest excited vibra- 
tional levels of the bare background-scattering potential to 
be -h x 0.0086 MHz and -h x 110 MHz, while the effec- 
tive range and resonance shift amount to r e s = 255 «o 
and Bq - fi res = -87 G, respectively. Accordingly, the 
model parameters of Eqs. (|A"2"9i l, (lA30l . (|A3TT >. and (TA321 
for the high-field zero-energy resonance are given by: i;\ = 
-6.5718 x 47T 3/2 fi 2 o- I /m, = -13.7782 x 4-7t 3/2 h 2 cr 2 /m, 
£ = 0.929532Ux2jT V4 h 2 /(mo- U2 \ <r x = 80.6241 a , cr 2 = 
38.512a , ando- = 18.7179 a . 
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